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Abstract. A large number of X-ray sources in globular clusters of galaxies other than the Milky 
Way has been found with Chandra. We discuss three issues relating to these sources. The X-ray 
luminosity function (XLF) of the sources in globular clusters of M31 is marginally compatible with 
the XLF of globular clusters of the Milky Way. The individual XLFs of a dozen elliptical galaxies, 
after correction for incompleteness, are compatible with one another and show no break; however, 
the XLF found by adding the individual XLFs of elliptical galaxies has a break at Lj; ~ 5 x 10^^ 
erg s^'. For the moment there is no evidence for a difference between the XLFs of sources inside 
and outside globular clusters of elliptical galaxies. It is not (yet?) possible to decide which fraction of 
low-mass X-ray binaries in elliptical galaxies outside globular clusters have formed inside globular 
clusters. 
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INTRODUCTION 

Observations with Chandra of an increasing number of nearby galaxies are revealing a 
sizable population of bright X-ray sources, many of which are in globular clusters. In 
elliptical galaxies, as many as half of all bright X-ray sources are in globular clusters; 
in spiral galaxies a smaller fraction of the bright sources is in globular clusters. The 
absence of recent star formation in elliptical galaxies implies that the X-ray sources in 
them are mainly low-mass X-ray binaries. Some questions that have arisen from early 
research are whether globular clusters do contain black holes? and whether the sources 
in elliptical galaxies outside globular clusters originate in globular clusters? 

Elliptical galaxies are a good target for the study of X-ray sources in globular clusters, 
because many have 5000 to 10000 globular clusters. The central galaxy of the Virgo 
cluster, M87, has 13500 globular clusters! The X-ray sources detected so far in the 
faraway galaxies are of necessity very bright, typically L^>5x 103^ ergs k Early 
studies indicate that about 4% of the globular clusters contains an X-ray source above 
this limit, roughly corresponding to 1 source per 5 x lO^L© / (Sarazin et al. 2003, Kundu 
et al. 2003). 

We have given an overview of the early papers on this topic in Verbunt & Eewin 
(2005). Comparison between studies is difficult, as they have different detection limits 
and compute X-ray luminosities in different energy bands based on different assumed 
spectra. The study of Kim & Eabbiano (2004) remedies this, and provides a first consis¬ 
tent comparison between a fair number (14) of different elliptical galaxies; we discuss 



FIGURE 1. Left: Comparison of the normalized cumulative X-ray luminosity functions (at Lx > 
ergs^^) of globular clusters of our Milky Way and ofM3L The Chandra luminosities given by Di Stefano 
et al. (2002) were multiplied by 0.46 to convert them to the energy range ofROSAT data from Verbunt et al. 
(1995). After Verbunt & Lewin (2005). Right: Comparison of the normalized cumulative X-ray luminosity 
functions (at Lx > ergs^^) of M87 sources in globular clusters and outside them. After Jordan et 

al. (2004). The Kolmogorov-Smirnovprobability that the normalized distributions are the same is 0.03 for 
the Milky Way vs. M 31, and 0.09 for in vs. out globular clusters ofM87. 


this important study in Sect. 3. In Sect. 2 we compare the X-ray luminosities of globular 
clusters in our galaxy and in the Andromea nebula. In Sect. 4 we discuss the question 
whether sources outside globular clusters can have formed inside them. 


GLOBULAR CLUSTERS IN THE MILKY WAY AND IN M31: 
DIEEERENT OR THE SAME? 

The Andromeda Nebula, a.k.a. M 31, was found with Einstein to have rather more bright 
X-ray sources in globular clusters than our Milky Way (Van Speybroeck et al. 1979). 
Since this discovery it has been debated whether this simply reflects the larger number 
of globular clusters of M 31, or an X-ray luminosity function with a higher fraction of 
bright sources in M31 (e.g. Verbunt et al. 1984). The debate continued in the ROSAT 
and Chandra era (Magnier et al. 1992, Supper et al. 1997, Di Stefano et al. 2002). In 
Figure[T]we show the cumulative X-ray luminosity functions for globular clusters of our 
Milky Way and of M 31. A 2-sided Kolmogorov-Smirnov test gives a probability of 3% 
that both distributions are the same. Thus, the evidence that the luminosity functions are 
intrinsically different is at the 2-sigma level: suggestive, but not conclusive. 

The X-ray luminosity functions of the globular clusters in M 31, and the bulge in M 31 
(Kong et al. 2002, 2003) are both roughly compatible with the X-ray luminosity function 
of elliptical galaxies (see next Section), an indication that the X-ray sources in all these 
old stellar populations are similar. (With a population of eight persistent sources and 5 
transients, the globular cluster system of our galaxy has too few bright X-ray sources to 
















constrain the slope of luminosity funetion well.) 


X-RAY LUMINOSITY FUNCTIONS OF GLOBULAR CLUSTER 
SYSTEMS OF ELLIPTICAL GALAXIES 

Kim & Fabbiano (2003) investigated the ineompleteness at the low-luminosity end of 
the X-ray luminosity function of the elliptical galaxy NGC 1316, and showed that the 
apparent break in the XLF disappears when appropriate eorreetions are made. They then 
investigated XLFs of 13 more elliptieal galaxies, and showed that none of these shows 
a significant break after correction (Kim & Fabbiano 2004, see also Gilfanov 2004). 
The elliptical galaxies comprise 7 galaxies of the Virgo eluster (at 17 Mpe), 2 of the 
Fornax eluster (19.9 Mpe), and 5 others (between 11 and 29 Mpe). The effeets leading 
to ineomplete deteetion effieieney at low X-ray luminosities are 

• the presence of diffuse emission of the hot interstellar medium in E and SO galaxies 

• the Eddington bias, enhaneing the number of sourees near the deteetion threshold 

• souree eonfusion 

• the larger point-spread function near the detector edge 

The X-ray luminosities are determined in the energy range 0.3-8.0keV; and sourees 
are eounted outside the innermost region of eaeh galaxy (r > 20") but within the 25 mag 
isophote. Correetions for the souree numbers and total X-ray luminosity of the galaxy 
are made on the assumption that they seale as the optieal luminosity, i.e. with this 
affeets the normalization, but not the slope of the luminosity funetion. It should be noted 
that Kim & Eabbiano (2003, 2004) do not diseriminate between sourees within and 
outside globular clusters, and give the XEEs of all sourees related to elliptieal galaxies. 
They find that the XEE of every individual elliptieal galaxy is compatible with a power 
law N{L^ oc with /3 = 2.0T0.2. The variation between the /3-values of different 
ellipticals can be aseribed to small-number statisties. In studies of individual elliptieal 
galaxies it is found that the XEE for sourees in globular elusters is the same as the XEE 
for sourees outside globular elusters (Maeearone et al. 2003, Sarazin et al. 2003, Jordan 
et al. 2004, see Eig.[T]). Thus, we may take the overall XEE as a proxy for the XEE of the 
globular eluster sourees. 

Sinee the individual X-ray luminosity functions appear to be the same, they can be 
added. The added XEE, at luminosities in the 0.3-8.0 keV band Lx> 10^^ ergs^\ 
of all 14 elliptieal galaxies is marginally eompatible with a single power law with 
/3 = 2.1±0.1.A broken power law gives a somewhat better fit, and is shown in EigureEl 
The best-fit break luminosity is (5.6±1.6)x 10^^ ergs“\ and the values for /3 at lower 
and higher luminosities are 1.8±0.2 and 2.8±0.6, respeetively. Remarkably, the break is 
at about the flux where the luminosity funetion of AGNs and quasars also shows a break. 
Whereas baekground sourees typieally eomprise about 5% of the sourees within the 
25 mag isophote, they are unlikely to eause the break in the added XEE of the elliptieal 
galaxies, whieh has been eorreeted for the baekground sourees (Kim & Eabbiano 2004). 

The break is at a luminosity whieh is about twiee the Eddington limit for hydrogen- 
rieh material, and eomparable to the Eddington limit for hydrogen-poor material (see 



FIGURE 2. Left: total corrected number of sources in 14 elliptical galaxies as a function of X-ray 
luminosity in the 0.3-8.0keV range. Right: corrected total X-ray luminosity of elliptical galaxies as 
a function of their absolute K-band magnitude Mk- The total X-ray luminosity is computed from the 
corrected X-ray luminosity function, extrapolated downwards to Lx = 10^^ erg The size of the symbols 

o scales with the number of globular clusters per unit luminosity of the galaxy • indicates a galaxy 
for which Sn is not known. (After Kim & Fabbiano 2004.) 


also Kuulkers et al. 2003). From their high X-ray luminosity and soft X-ray spectrum, it 
has been argued that the sources above the break are accreting black holes (e.g. Angelini 
et al. 2001). As regards the sources below the break, Bildsten & Deloye (2004) note 
that the assumption that most very luminous X-ray sources are ultracompact binaries, in 
which the donor star is a hydrogen-depleted degenerate star, explains the location of the 
break at the Eddington limit of hydrogen-poor material, the high X-ray luminosities, the 
source incidence (or birth rate), and the X-ray luminosity function. The model does not 
explain why sources outside globular clusters, which are less likely to be ultracompact 
binaries (if our own Milky Way is a guide), have the same X-ray luminosity function. 
This brings us to the question whether sources outside globular clusters could have 
formed inside them. 


OUTSIDE AND INSIDE GLOBULAR CLUSTERS 

In elliptical galaxies it is found that the spatial distribution of X-ray sources outside 
globular clusters is similar to the spatial distribution of globular cluster X-ray sources 
(e.g. Jordan et al. 2004 for M 87). Similarly, the X-ray luminosity functions of globular 
cluster sources and other sources are very similar (Fig.[T]). This has led to the revival of a 
suggestion by Grindlay & Hertz (1985) that all low-mass X-ray binaries, including those 
now outside globular clusters, were formed inside globular clusters (White et al. 2002). 

In the Milky Way and in M 31 there are about 10 bright low-mass X-ray binaries in 
the disk for each one in a globular cluster. In elliptical galaxies, there is about 1 source 
outside globular clusters for each one in them. This suggests that the majority of low- 

















mass X-ray binaries in the disk of the Milky Way , M 31 and by generalization in spiral 
galaxies are formed in the disk, and not in globular elusters (Verbunt & Lewin 2005). It 
should be noted that we eannot compare fractions inside and outside globular clusters 
for the Milky Way and M 31 on one hand and elliptical galaxies at the other hand at the 
same luminosities, because there are very few X-ray sources in the Milky Way and M 31 
above ~ 5 x 10^^ erg s~\ which is the lower limit of detectable sources in ellipticals 
(see Fig.[T]). 

The total X-ray luminosity of low-mass X-ray sources in elliptical galaxies scales with 
the luminosity (hence presumably stellar mass) of the elliptical galaxy (Figure|2l). This is 
expected for sources formed outside globular clusters. It is also expected if the sources 
are mainly formed in globular clusters, because the number of globular clusters is on 
average higher in large elliptical galaxies than in small ones. Correlations of the total X- 
ray luminosity have been made with the specific frequency Sm of globular clusters (i.e. 
the total number of globular clusters divided by the luminosity of the galaxy). We refer 
to this here as the "global" specific frequency which includes all globular clusters in the 
entire galaxy. These correlations greatly suffer from a serious problem related to our 
lack of knowledge of the global Sm- Reliable values for the number of globular clusters, 
in general, come from HST-WFPC2 with a very small (5.7 square arc minutes) field of 
view. This provides only a reliable local value for Sm but not the global value (see the 
discussion in Section 8.3 of Verbunt & Lewin 2005). 

Kim & Fabbiano (2004) note that the correlation between the total X-ray luminosity 
(of low-mass X-ray sources) and the infrared luminosity of elliptical galaxies has a scat¬ 
ter which is larger than the measurement accuracy. In Fig. El we indicate the elliptical 
galaxies with a large St^ with a large symbol. We see that the elliptical galaxies brighter 
in X-rays than the average correlation indeed have a tendency to have a higher This 
indicates that globular cluster sources contribute significantly to the total X-ray luminos¬ 
ity of the galaxy; to argue that sources outside globular clusters originate inside globular 
clusters, one would have to show that the total X-ray luminosity (or better, number) of 
the sources outside globular clusters separately scales with specific frequency Syy. In the 
sample of galaxies studied by Kim & Fabbiano (2004) there is no correlation between 
Mjc and (see Fig.EJ). 

CONCLUSIONS AND OUTLOOK 

The conclusions drawn above, i.e. that the X-ray luminosity function of the X-ray 
sources in old populations is universal and has a break dXLx — 5x 10^^ erg s^^ are based 
on a study which does not discriminate between sources inside and outside globular 
clusters. With the increasingly large sample of sources it is possible and necessary 
to repeat this study for the globular cluster sources separately. This may also help 
in adressing the question which fraction of the sources outside globular clusters were 
actually formed inside them. The location of the break suggests that the sources above 
the break are accreting black holes. The sources below the break may be dominated by 
ultracompact sources (Bildsten & Deloye 2004). 

There are two important areas of further study not discussed above. The first is the 
question why clusters with a high metallicity have a higher probability of containing 



FIGURE 3. The number of neutron stars per solar mass formed from an initial mass distribution 
N{m) as a function of the index a. The initial range of stellar masses is taken to be from 0.08 

to 80 Mq, and of neutron-star progenitors from 8 to 80 Mq. The solid curve is normalized on the initial 
mass of the cluster, the dashed curve on the current luminous mass, assuming that stars from 0.8 to 8 Mq 
have evolved into white dwarfs. 


an X-ray source, even if the number of eollisions in them is the same. An interesting 
ealeulation by Jordan et al. (2004) shows that a moderate dependenee on metallieity of 
the slope of the Initial Mass Function leads to a differenee in the numbers of neutron stars 
in metal-rieh and metal-poor clusters (per unit mass) which is large enough to explain 
the observed preferenee for high-metallieity elusters (see Fig.|3l). Briefly, if we write the 
number of stars at mass m =M/Mq as N{m) = Km^^, we have for the number of stars 
with initial mass between mi and m 2 (that have evolved into neutron stars) per unit mass 
between nia and nit,: 


N{nil,m 2 ) 2 —a m} “ —m^ “ 

M{ma,mb){MQ) \-a 

Jordan et al. (2004) discuss the number of neutron stas per unit of initial mass of the 
eluster, i.e. they take mt, = m 2 . This ratio, for values m^ = 0.08, m/, = m 2 = 80 and 
mi = 8, is shown as a solid line in Fig. E] as a funetion of a. If a is lower at higher 
metallieity, i.e. metal-rich clusters have a shallower initial mass funetion, then metal- 
rich clusters have a relatively higher number of neutron stars. We may add that this 
effeet is even stronger if we seale not on the initial mass, but on the eurrent luminous 
mass, i.e. the stars that determine the observed K magnitude. Eq.[T] for m 2 = 80 as 
before, but with = 0.8, is shown in Fig.|3l as a dashed line. (Stars in the range 
0.8 < m < 8 have evolved into white dwarfs.) This deserves further study, in partieular 
whether the required dependenee of the slope of the IMF on metallieity ean be more 
firmly established. 

The seeond researeh area is the indieation that the probability for a globular eluster 
to contain an X-ray source increases slower with density than the number of eollisions 






occurring in the cluster (Jordan et al. 2004). The models for formation and evolution of 
X-ray sources in globular clusters must address these questions. 
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